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^ ' Abstract 

. We study the Landau gauge photon propagator in the three-dimensional AbeUan Higgs 

, model with compact gauge field and fundamentally charged matter in the London limit. 
The total gauge field is split into singular and regular parts. On the confinement side of 
the string breaking crossover the momentum dependence of the total propagator is char- 
, acterized by an anomalous dimension similarly to 3-D compact QED. At the crossover and 
throughout the Higgs region the anomalous dimension disappears. This result perfectly 
agrees with recent observations that the monopole-antimonopole plasma leads to nonzero 
anomalous dimension and the presence of the matter fields causes monopole pairing into 
I dipole bound states. The Yukawa mass characterizing the propagator part from regular 
gauge fields is non-vanishing at the Higgs side and coincides with the mass found for the 
total propagator. The regular gauge field without anomalous dimension becomes massless 
at the crossover and in the confinement region. 



1 Introduction 



Nowadays, the interest in the lattice Abehan Higgs model with compact gauge field (cAHM) 
in three dimensions has grown because of its relation to high energy physics |jl|, and its 
apphcations in condensed matter physics [Q. 

The compactness of the gauge field leads to the presence of monopoles which are instanton- 
like excitations in three space-time dimensions. Being in the plasma state the monopoles 
and antimonopoles of this theory guarantee linear confinement of electrically charged test par- 
ticles [Q. They are forming an oppositely charged double sheet along the minimal surface 
spanned by the Wilson loop (i.e. the trajectories of the heavy charges). Due to screening, the 
free energy of the surface increases proportionally to the area of the surface such that an area 
law for the Wilson loop emerges. 

Within the confinement region of cAHM corresponding to small couplings of matter fields to 
gauge fields the monopole-antimonopole plasma state is still realized. As the hopping parameter 
increases, the system enters the Higgs region where monopoles and antimonopoles become 
bound into magnetically neutral dipoles. This scenario has been demonstrated in a preceding 
paper |^ and it has been related to the phenomenon of string breaking. Indeed, when monopole 
pair formation occurs, this results in the breakdown of linear confinement at large distances. 
Usually this is interpreted in an alternative way saying that dynamical matter fields in the 
same representation as the external test charges break the confining string by screening the 
charge of the latter. This argument is applied, irrespective whether the dynamical matter field 
is fermionic (the quarks in QCD) or bosonic (the Higgs particles in our case). Let us stress here 
the other point of view, according to which the monopole mechanism of confinement is changed 
in a way to produce a different form of the inter-particle potential. At large separations 
R of the charges the string tension should be absent. However, if i? -C Rhv (where i?br is 
the characteristic string breaking distance) the test charges are able to recognize individual 
monopoles even if they are bound in dipoles. Therefore, the monopole and antimonopole fields 
may induce a piecewise linearly rising potential. These simple considerations can be made more 
rigorous by analytical calculations for a gas of dipoles with small magnetic moments. 

In order to destroy the linearly rising potential within CAHM3, the coupling between Higgs 
and gauge fields must be sufficiently strong. It would be tempting to associate the onset of 
string breaking with a phase transition between confinement and Higgs phases. However, it 
was numerically shown that in the London limit of CAHM3 the string breaking happens in 
a region of the phase diagram where a first or second order phase transition can definitely be 
excluded. In the present paper we will call this the "string breaking crossover" , but one should 
keep in mind that a thorough reconstruction of the monopole configurations is accompanying 
this. However, the possibility [§] that string breaking (and monopole pairing) is associated with 
a Berezinsky-Kosterlitz-Thouless (BKT) type transition is not ruled out. 

Recently, it was found that the matter fields in the Abelian Higgs model lead to a logarithmic 
attraction between monopoles and antimonopoles which would explain the formation of 
monopole-antimonopole bound states and string breaking. Adding massless quarks also forces 
the Abelian monopoles and antimonopoles to form bound states [§. Note that the origin of 
monopole binding in the zero temperature case of the cAHMs is physically different from the 
monopole binding observed at the finite temperature phase transition in compact (2 + 1)- 
dimensional pure QED ||^, |1^ . 
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We have recently studied the effect of finite temperature deconfinement of (2+l)-dimensional 
cQED on the photon propagator in Refs. |[ll|, |l^.We could demonstrate that the momentum 
behaviour of the photon propagator in this theory is described, rather similar to gluodynamics, 
by a Debye mass and by an anomalous dimension which both vanish at the deconfinement 
transition. This mechanism could be clearly attributed to pairing of magnetic monopoles. The 
monopole-antimonopole plasma contribution is relatively easy to exhibit by explicit calculation 
and can be eliminated by monopole suhtmctioniiom. the ensemble of gauge fields. 

As for gluodynamics, which motivated our study of the CAHM3, numerical lattice results 
show that the propagator for all these gauges in momentum space is less singular than l/p^ 
in the immediate vicinity of = 0. Recent investigations in the Landau gauge and in the 
Laplacian gaugeQ show that, beside the suppression at — > 0, the propagator is enhanced 
at intermediate momenta which can be characterized by an anomalous dimension |15|. This 
enhancement of the Landau gauge propagator in SU{2) gluodynamics has been interpreted ||T6[ 
by focusing on P-vortices appearing in the maximal center gauge. Subtracting the vortices 
removes the enhancement at intermediate momenta. The results for the propagator at zero 
momentum are ranging from a finite |[T^ (Laplacian gauge) to a strictly vanishing |T^, |TS| 
(Coulomb gauge) value. The vanishing of the Landau gauge propagator at = 0, suggested 
by considering the Faddeev-Popov mechanism, remains obscured so far in the results, probably 
because the lattices are still too small. 

In this paper we are going to investigate in which way the photon propagator within cAHMa 
changes at the string breaking crossover, turning from confinement to the symmetry-broken 
Higgs region. As in our previous work we have chosen the propagator in the minimal Landau 
gauge because this is the covariant gauge which has been adopted in most of the investigations 
of the gauge boson propagators in QCD [l^, |l^] and QED |jl9|. For the behaviour at the 
string breaking crossover of CAHM3 we anticipate that a confining propagator will change into 
a Yukawa-like propagator corresponding to the onset of the Higgs mechanism. 

The present paper is structured as follows. Next, in Section ^, we will recall the model 
and the definition of the photon propagator. The form of the fitting function and the method 
of monopole subtraction are also introduced there. In Section ^ we will report the numerical 
results of the present study and discuss the sensitivity with respect to the Gribov copy problem. 
We conclude in Section ||. 



2 The Model and the Definition of the Propagator 

We consider the three-dimensional Abelian gauge model with compact gauge fields d^^fj. living 
on links and a fundamentally charged Higgs fields on sites. For simplicity we consider the 
London limit of the model, which corresponds to an infinitely deep potential for the Higgs field. 
Consequently, the radial part of the Higgs field, \^x\, gets frozen. With (p^ = arg'^'x we define 
the model by the action 

5* [61] = -pJ2cos9p - kJ2 cos{9x,^, + (fx+fi - ^x) , (1) 

P 

^In order to avoid the problem of Gribov copies the alternative Laplacian gauge has also been used 
recently jlj] in Yang-Mills theory. However, principal questions of renormalizability and transversality are yet 
unsolved in this case. 
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where (3 is the inverse gauge couphng squared, k is the hopping parameter, and 6'p is the 
plaquette angle representing the curl of the link field 6x,^i- 



For the simulations we use a Monte Carlo algorithm similar to the one described in Refs. |^ 
12] , ^. The Higgs field angles have been updated in alternating order with the gauge field angles. 
In both cases, one 5-hit Metropolis sweep together with 2 microcanonical sweeps constitute a 
total gauge or Higgs update. Global updates of the gauge field have a negligible acceptance rate 
for K 7^ 0. They have been discarded at all. The numerical calculations have been carried out 
for fixed gauge coupling, [3 = 2.0, on lattices of size 32^, as in Ref. 0. For these parameters, 
the low-K region ("confinement region") and the high-K region ("Higgs region") of the phase 
diagram are separated by a crossover. The crossover point — observed as a position of the 
peak of the susceptibility of the hopping term Sgh = J2x,^l cos{9x^^ + fx+fi — V^x) — is located 
at Kc = 0.526(1). We have considered from 0(200) to 0(400) independent configurations 
(obtained after 10 subsequent updates) to measure the propagator. 

The discussion of the photon propagator (and its various parts) is given in lattice momentum 
space. Being always defined in the context of a specified gauge, in our case the minimal Landau 
gauge 

Ecos(0^max (2) 

X,fJ, 

with respect to gauge transformations G, the propagator is written in terms of the Fourier 
transformed gauge potential, 

which is a sum over a set of points x = n + ^fl, the midpoints of the links in /i direction, 
which form the support of A^^^ on the lattice, n denotes the lattice sites (nodes) with integer 
Cartesian coordinates. The propagator is the gauge-fixed ensemble average of the following 
bilinear in A, 

D,M = (hAO ' (4) 

where the lattice momenta p on the left hand side of (^ are related to the integer valued Fourier 
momenta k by the expression (a is the lattice spacing): 

p,(M = -sin^, A;, = 0,±1,...,±:^. (5) 

The lattice equivalent of = is in 3 dimensions 

v\k) = ^^jzU^'^ . (6) 

For this paper, we decided to identify the gauge field A$^y, in terms of the sine function of 
the link angle (with g\ = l/(a/5)) 

A^^^ = -^smOn,^. (7) 
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The corresponding propagator, called DJ^^ in Ref. [|T2|, will be called simply D^i, for brevity. The 
extraction of the Fourier transformed gauge field (and of its regular and singular components) 
has to be performed after the original gauge field configuration has been put into the minimal 
Landau gauge. The procedures employed for gauge fixing in the present context have been 
described at length in Ref. . 



At zero temperature, for perfect Euclidean rotational invariance, the continuum propagator 
would be expressible by functions of p^. The most general tensor structure is then the following 
one including two scalar functions of p^, 

D,Ap} = P,Ap}D{p^) + P^^ (8) 

with the three-dimensional transverse projection operator 

P,uip) = V - ^ • (9) 

The two structure functions D{p'^) and F{p^) can be extracted by projection, on the lattice 
from D^y{p) according to (^, as 

3 

F{P^) = E D,,{p) (10) 

and 

p'D{p') = l j: P,Ap)D,Ap}. (11) 
^,1^=1 

They are approximately rotationally invariant, i.e. individual momenta p might slightly differ 
in the function values D or F they provide, even if they have the same p"^. Dense, in p"^ nearby 
data points may scatter rather than be forming a smooth function of p^. 

In practice, using these definitions, we extract at first the function F{p'^) from Eq. ([To|). In 
the sum the imaginary parts of non-diagonal D^^ cancel. Then, D{p'^) is obtained through 



D{p') = - D,,{p) + D,,{p) + D,,{p) - . (12) 




For exactly fulfilled Landau gauge F{p'^) = 0. On the lattice, in the case of the sine- 
definition used for A^^^ [Eq. (|^)], this is actually the case as soon as one of the local maxima of 
(0) (Gribov copies) is reached, with an accuracy which directly reflects the precision at stopping 
of the gauge fixing iterations. 

The effect of monopoles (singular fields) can be distinguished from that of the regular 
( "photon" ) fields using the splitting of the gauge field angles d^^fj. into a regular and a monopole 
part following Refs. |]20| , |Tl| , |l^. In the notation of lattice differential forms this can be written 
as: 

e = 6'^^ + ^"^"^ , e""""" = 27iA-^6p[j] , (13) 
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where is the inverse lattice Laplacian and the 0-form *j G ^ is nonvanishing on the sites 
of the dual lattice occupied by monopoles and antimonopoles. The 1-form *p[j] corresponds 
to Dirac strings (living on the links of the dual lattice) which connect monopoles with anti- 
monopoles, = *j- For any Monte Carlo configuration, we have fixed the gauge, then 
located the Dirac strings, p[j] ^ 0, and constructed the monopole part 6*™°"^ of the gauge field 
according to the last equation in (|l^). The regular photon field^is taken just as the complement 
to the monopole part according to the first equation of (|13]). 

The regular and the singular parts of the gauge field contribute to the propagator as follows: 
the total propagator decomposes like D = D^^^ + D^°^ + D'^^^, where D"^^^ represents the mixed 
contribution from regular and singular fields. An analogous decomposition is valid for the 
longitudinal structure function F which vanishes to a good accuracy for the total propagator. 
In this paper we are interested mainly in the regular part and in the total photon propagator. To 
describe the propagators quantitatively we have fitted the total propagator using the following 
function: 

D(p') = jj^msr^;^^)+c, (14) 

where Z, a, m and C are the fitting parameters. This fit has been successfully used to describe 
the propagators of the finite and zero-temperature compact U(l) gauge model [0, |12[ in 2 + 1 



or 3 dimensions, respectively. The form is similar to some of Refs. [|T5l where the propagator 



in gluodynamics has been studied. The meaning of the fitting parameters in Eq. (JT^ is as 
follows: Z is the renormalization of the photon wavefunction, a is the anomalous dimension, m 
is a mass parameter. As shown in Ref. [|12], in CQED3 this mass parameter coincides with the 
Polyakov prediction ^ for the Debye mass, generated by the monopole-antimonopole plasma. 
The parameter C corresponds to a 5-like interaction in the coordinate space and, consequently, 
is irrelevant for long-range physics. 

The regular (or "photon") part of the propagator, D'^^^, has been fitted by the Yukawa form 
(with an additional contact term): 

/9 (p2 + m2^g) 

The fits together with the appropriate propagator data will be presented in the next section. 



3 Numerical Results 

We begin with the typical shape of the total photon propagator and its regular part, respectively, 
at fixed (3 = 2.0. The total photon propagator D, multiplied by p^, is shown in Figure |^(a) for 
various values of k as a function of We also present the fits of the data by the function (|I^) . 
One clearly observes that the total propagator is well described by the fitting function in 
each particular case. The "infrared" enhancement at intermediate momenta disappears with 

^In principle, the regular part of the links could have been reconstructed as well, without recurrence to 
the singular part. In this case the regular propagator would become completely independent of the number of 
Gribov copies Nq (see the discussion below). 
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Figure 1: Momentum dependence of the total propagator (a) and its regular part (b) and 
fitted curves 



increasing k. In contrast, the suppression at very small |p| remains at all n values. We can 
conclude that the total photon propagator is less singular in the infrared than a free one. 
According to our fitting function ([T^) the propagator is finite at = 0. 

For comparison, the regular photon propagator part D^'^^, multiplied by p^, is shown in 
Figure |I](b) for various values of k. One can see that the behaviour of D^'^^ oc is only 
observed in the confinement region; in the Higgs region the propagator becomes massive. The 
fits using function (|T5|) work very well everywhere except for smallest non-zero momenta in the 
closest vicinity of the crossover point. 

Note that the data for both the total propagator and the regular part have been averaged 
over lattice momenta corresponding to the same before fitting (as in Refs. [^, 0). 

From the fits we have obtained the characterizing fit parameters for the transverse photon 
propagator D in gauge fixed ensembles with various numbers {Nq < 100) of Gribov copies. 
This means that, after applying the gauge fixing algorithm to Nq random gauge copies of 
the original Monte Carlo gauge field, only the configuration with the maximal value of the 
gauge fixing functional (H) is kept for evaluation. In Fig. ^ we show the dependence of the 
fit parameters for the total photon propagator on the number of Gribov gauge copies. Three 
typical K values are used: k = 0.3 in the confinement region, k = 0.535 very near the crossover 
and K = 0.7 in the Higgs region. While a moderate number Nq ~ 30 seems to be sufficient for 
convergence of all parameters at all k values, there is no clear tendency between different k's. 
The parameters describing the regular propagator D'^'^^ converge within the first few {No ~ 5) 
copies. 

To keep the influence of Nq on the propagator negligible we have used in the final measure- 
ments Ng = 60. The resulting fit parameters of interest are presented in Fig. ^ as functions 
of K. Fig. |^(a) depicts the mass taken from the fit (|15[) for the regular part of the propaga- 
tor ("regular mass" nij-eg), and from the fit (p^ for the total propagator ("total mass", m). 
The dashed line represents the expected behaviour of the mass in the limit of large k and (3 
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Figure 2: Dependence of the fit parameters m, a, Z and C on the number Nq of Gribov copies, 
couphngs: 

mth(/3, 



(16) 



where Pv and ny are the Villain couplings corresponding to the Wilson couplings /? and k. 
This parameterization is analogous to Ref. PT|, with 



Pvm= 2 log 



-1 



2 log 



(17) 



Note, that although (3 = 2.0 (used in our calculations) is not very large, the first relation in 
Eq. ([T7|) works with an accuracy of a few percent already at /? ~ 1 according to the calcula- 
tion |2l| of the critical coupling in AD cQED. 

We observe that both the total and regular masses coincide with each other on the Higgs 
side of the crossover, and there they are very close to the prediction (|16|,|l^, the immediate 
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vicinity of the crossover point excluded. This observation can be easily understood. Two 
different sources contribute to the gauge boson mass: one is arising nonperturbatively from 
the monopoles (the Debye screening mass) while the other is due to the explicit presence of 
the mass term in the action (|lD of the model. The Debye mass generation works in the case 
of the monopole-antimonopole plasma and it is obviously absent in the magnetic dipole 
gas . Since at large k monopoles and antimonopoles are bound into dipoles ^ , the Debye 
contribution to the mass disappears, and the gauge boson mass is exclusively given by Eq. (|16|). 
As can be seen from Fig. 0(a), at small k values the total mass m is close to the Debye mass 
of pure CQED3 {i.e. to the mass value at k = 0) because the effects of the mass term in the 
action (^ are small. Moreover, in accordance with our expectations, the regular mass rrireg 
vanishes. 
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In the region very close to the crossover, k 



the mass m shows a minimum caused, 



as one could guess, by the interference between perturbative mass and Debye mass effects. 
Indeed, as k, tends to Kc, the Debye mass gets smaller since the density of the monopole- 
antimonopole plasma drops rapidly. One the other hand, the perturbative mass term becomes 
more significant. The interplay of these two tendencies results in the noticed minimum at 

In Fig. ^(b) we show the anomalous dimension a for the total propagatoiQ as a function of 
the hopping parameter n. The anomalous dimension is non-vanishing in the confined region 
and it turns to trivial values at the crossover]^, a ^ 0. This behaviour can be compared with 
our studies of compact QED2+1 at finite temperature ||lT|, |12[: in the monopole-antimonopole 
plasma phase (corresponding to the confinement side) a > while in the magnetic dipole phase 
(corresponding to the Higgs side) a = 0. Thus, in cAHM the effect of the monopole pairing 
on the propagator is the same as in cQED: the anomalous dimension gets close to zero in the 
Higgs phase dominated by the magnetic dipole gas. 

Similarly to the anomalous dimension, the effect of the monopole pairing on the renormal- 
ization parameter Z of the photon wavefunction in cAHM, shown in Figure |^(c), is remarkably 
similar to the cQED case observed^ in Refs. |TT| , |T^. The total photon factor Z suddenly drops 
at the crossover point, n = Hc while the regular factor, Zj-eg, is almost insensitive with respect 
to the crossover {Z Z^-eg ~ 1 for all k values). The small contact term parameter of the total 
propagator, C, changes its sign at the crossover. Figure ||(d). The small and negative contribu- 
tion of Creg in the confinement is responsible for the tiny decrease of the massless p^D'^^^ with 

[seen at zero and small k in Fig. |^(b)] . A zoom to the crossover region is shown in Fig. ^ 




Figure 4: Zoom to masses and anomalous dimension in the crossover region. 

^Recall, that the anomalous dimension for the regular part of the propagator is zero, Eq. ([l5|). 

*The 4-parameter fits lead even to small negative a in the region k 1. However, fitting there the propagators 
with fixed a = 0, the obtained masses coincide within errors with those of the unconstrained fits. So we associate 
this behaviour with statistical fluctuations. 

^ Due to an error, the Z factors in Figs. 5c and 7c of Ref. [|2| need to be corrected by a factor to" before 
they can be compared with the Z given here at k = 0. 
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where all effects described above for m, mreg and a are seen more clearly together. 



4 Conclusions 

In our numerical studies we found that the gauge boson propagator in the London limit of the 
three-dimensional compact Abelian Higgs model possesses a non-zero anomalous dimension 
below the string breaking crossover. The effects of the matter fields on the propagator are 
remarkably similar to the finite-temperature deconfining effects in the pure gauge compact 
[/(I) model. Both in the Higgs (deconfinement) region of cAHM and in the deconfinement 
(high temperature) phase of cQED the anomalous dimension of the propagator vanishes. In 
the confining regions of both models the anomalous dimension is non-vanishing. The positive 
anomalous dimension is due to presence of the monopole-antimonopole gas in the plasma 
state |TT|. As we move towards the deconfinement phase a decreases and becomes zero when 



the monopole-antimonopole plasma turns into a magnetic dipole gas. However, the origin of 
the pairing phenomenon in both models is different: in cAHM the monopole pairing is caused 
by the matter fields while in cQED the monopoles form the bound states due to temperature 
effects in the monopole-antimonopole interaction. 

For the limit of cAHM with frozen radial Higgs degrees of freedom we have found that the 
anomalous dimension for the gauge field does not become clearly negative together with the 
onset of string breaking. However, if radial fluctuations of the Higgs field would be allowed, 
the emergence of a negative anomalous dimension cannot be ruled out |^. 



The mass in the gauge boson propagator closely follows the tree level expectations on the 
Higgs side of the transition and simply corresponds to a massive Yukawa propagator practically 
unaffected by remaining dipoles in the vacuum. In the confinement region the mass related to 
the total gauge field is of Debye type while the mass of the regular or photon part of the gauge 
field degrees of freedom remains massless. 
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